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ABSTRACT 

We present an observational test of the spectrum of the ionizing background at ^ ~ 5; 
the test is sufficiently sensitive to determine whether Popii or Pop ill stars are the 
dominant source of ionizing radiation. The ionizing background at 2 ~ 5 may reflect 
the nature of the sources responsible for the final overlap phase of reionization. We 
find that rest-frame extreme-UV Hei absorption will be detectable in deep spectral 
observations of some rare z ~ 5 quasars; the ratio of He i to Hi absorption reflects 
the shape of the ionizing background in the photon energy range between 13.6 and 
24.6 eV. Most z ~ 5 quasars have too much Hi absorption along their line of sight 
for He I absorption to be observed. However, based on current measurements of Hi 
absorber statistics, we use Monte Carlo simulations to demonstrate that the Sloan 
Digital Sky Survey (SDSS) will discover a sufficient number of 2; ~ 5 quasars to 
turn up a quasar suitable for measuring Hei absorption (and we illustrate a selection 
method to identify that quasar). From simulated observations of a suitable 2 ~ 5 
quasar with a 10-meter telescope, we show that a constraint on the spectral slope of 
the ionizing background at that redshift can be obtained. 

Key words: quasars: absorption lines - intergalactic medium - diffuse radiation - 
cosmology: theory 



1 INTRODUCTION 

InterRa lactic hydrogen is almost completel y io nized by 
z ~ 6 jBecker et alJl200ll : ^jorgovski ct al. 2003), but the 
sources re sponsibl e for rei onization ar c still not known (see 
iBarkana fc Loeb 200ll and lLoeb fc'^rkana 2001 for reviews 
of the theoretical possibilities). A low neutral- hydrogen 
fraction is found at lower redshift even in relatively over- 
dense systems, despite the short (compared to the Hub- 
ble time) recombination time-scale in those systems. This 
implies a continued supply of ionizing photons after reion- 
ization completed. Measurements o f the intensity of the 
ionizing background below a ~ 3 iScott et al.l I2OOOI) are 
consistent with an ionizing background produced by ob- 
served quasars jHaardt fc Madaul Il996l) . At higher red- 
shifts . however, the de clining abundance of bright quasars 
(e.g., iFan et alJl2001b^ suggests that they cannot provide 
an ionizing b ackground sufficient to reionize the universe 
iWvithe fc Loe b 2003 a), and this conclu sion is supported 
by deep x-ray surveys llBarger et alJl2003ll . 
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Recent results from the WMAP satellite imply that 
reionization was s ubstantially underw ay by z ~ 15 
jKogut et al.l l2003l : ISoergel et all l2003|). However, quasar 
absorption stud ies at 2 ~ 6 ^Becker et al.l I2OOII: 
iDiorgQvski et alJl200ll) suggest a dramatic evolution in the 
ionizing background around that redshift. Some authors 
have reconciled the WMAP evidence for early reionization 
with the quasar result that reionization is just finishing 
at 2 ~ 6 by postulating two reioniz ations of the universe 
iCenll2003alrbt IWvithe fc Loebll2003allbl) . generally invoking 
a hard Pop ill stellar spectrum for the early reionization and 
a Pop II stellar spectrum for the final reionization at 2 ~ 6. 

The most direct technique to ascertain the sources re- 
sponsible for the ionizing background at high redshift is to 
take a census of all sources that contribute ionizing photons. 
The advantage of this technique is that the relative contri- 
bution of different types of sources are measured directly, 
as are the spatial distributions and other properties of the 
populations. Unfortunately, the large luminosity distance to 
high redshift means that only the most luminous sources can 
be detected wi t h curre nt observational methods. For exam- 
ple, iFan et alJ ^2001b^ measure a power-law high-redshift 
quasar luminosity function with a steep slope, but estimate 
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that the traditional power-law break of the quasar luminos- 
ity function would occur more than one magnitude below 
their survey limit: there is little constraint on the faint end 
of the quasar luminosity function from direct source detec- 
tion. 

A complementary method to direct detection is mea- 
surement of the properties of the integrated ionizing back- 
ground. Although this technique does not identify ionizing 
sources directly, it is effective at detecting ionization due 
to a very abundant population of low-luminosity sources. 
Moreover, the shape of the integrated ionizing background 
is dictated by the combination of sources responsible for 
it; different spect ra arise from vario us possible contributors, 
such as quasars llTelfer et ahllioOSl) . Popii OB stars (e.g., 
iLeitherer et a l*199^. very massive metal-free (Pop ill) stars 
lEromm^ Kudritzki & Locb 2001), and x-rays from early 
stars f'Qhl l2'00lh . 

This paper describes possible observations that would 
constrain the spectral shape of the ionizing background just 
after reionization finished. In Section |5| we introduce quasar 
absorption line spectroscopy in the context of this paper. 
Section|3|describes our parametrization of the ionizing back- 
ground, and the resulting ionization state of absorption sys- 
tems in the intergalactic medium (IGM). In Section |1| we 
apply those results to observed statistics of Hi absorbers 
to generate a model for Hi and He I opacity toward high- 
redshift quasars; the assumed properties of the quasars are 
given in Section |S] The method of simulating line of sight 
(LOS) absorption spectra and studying the sensitivity of the 
results to the shape of the ionizing-background spectrum is 
described in Section |S] Section Q presents a method to se- 
lect the best quasars for absorption line spectroscopy. We 
summarize our results in Section |S| 
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Figure 1. Simulated absorption spectra for 2 = 5 quasars with 
ABi45o = 20. The top axis is labelled in rest-frame wavelength, 
and the bottom axis is labelled in the corresponding observed 
wavelength at 2 = 5. The top panel shows the wavelengths where 
the first four Lyman series line forests begin; it also shows where 
the the first two ground-state He I transition opacities begin. The 
He: band is shown, as is the region of the corresponding Lya 
forest. The middle panel shows a quasar spectrum processed by 
IGM absorption, along a line-of-sight with relatively little IGM 
opacity. The bottom panel shows a quasar absorption spectrum 
along a typical line-of-sight. 



2 ABSORPTION SPECTROSCOPY 

The easiest way to make a crude measurement of the shape 
of the ionizing background is through the ionization state of 
photoionized IGM gas; in the next Section we show explic- 
itly how the background spectrum is related to the ionization 
state. Atoms and ions with at least one bound electron in 
the ground state have strong bound-bound transition cross- 
sections, typically located at UV or soft x-ray energies. Con- 
sequently, even small amounts of these species on the line 
of sight to a background source produce strong absorption 
features. The most common extragalactic application of this 
technique uses UV-bright quasars as the background sources 
(for a recent review see IBechtold .200311 . Many species have 
been detected in absorption, ranging from Hi and Hell, to 
Fe II and Zn ll, plus high-ionization ions such as C iv and 
Si IV (e.g.,^chaska et al. 2001). 

We examine the first ionization states of hydrogen and 
helium; they are the most abundant elements in the IGM, 
and are responsible for most of the UV opacity toward high- 
redshift quasars. Specifically, we will investigate the relative 
neutral fractions of hydrogen and helium, which have first 
ionization potentials of 13.6 and 24.6 eV, respectively. Since 
the photoionization cross-sections for H I and He I are peaked 
at threshold, their relative neutral fractions constrain the 
effective slope of the ionizing background between 13.6 and 
24.6 eV. 



The relative abundances of Hi and Hel constrain the 
ionizing background — the challenge is to measure these 
abundances in the IGM at high redshift. Both atoms ex- 
hibit line and continuum absorption. For H I, the Lyman se- 
ries lines begin at 1216 A and continue down to 912 A, where 
continuous photoelectric absorption begins. The correspond- 
ing line transitions in He I are from Is^ to ls2p (584 A), ls3p 
(537 A), etc., down to the first ionization threshold of he- 
lium, at 504 A. These wavelengths are the rest-frame values; 
the observed wavelengths of these transitions depend on the 
redshift of the absorbing gas. If the neutral gas is distributed 
in many systems with a discrete redshift distribution, then 
the associated absorption lines appearing in a quasar spec- 
trum form a 'forest' of absorption features against the con- 
tinuum. Figure0shows two example absorption spectra with 
the general absorption regions labelled. 

The absorption cross-sections of H I and He I are larger 
at the center of a line (at the IGM temperature of ~ 
2 X 10'* K) than in the continuous absorption region. Despite 
this fact, the strongest absorption effect for the expected 
distribution of neutral IGM gas (see Sectional is from the 
cumulative continuous absorption due to many absorbers at 
different redshifts, named the 'valley' by M0llcr & Jakobse^ 
il990). This is illustrated for Hi in Fig.0 once Lyman limit 
absorption begins (shortward of 912(1 -I- zq) A, where zq is 
the quasar redshift), the forest typically gives way to al- 
most complete absorption. In such cases it is very difficult 
to observe helium forest lines because of the small residual 
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continuum flux at short wavelengths. In rare cases the valley 
is not very pronounced, and the continuum is strong enough 
that the He I 584 A forest region may be observed. We call 
the region of the He I 584 A forest red- ward of He I 537 Athe 
'He I band.' 

The identification of absorption features in spectra like 
those in Fig. can be difficult, especially in low-resolution 
spectra, except in the region red- ward of the Ly/3 forest, 
where the only possible absorption feature from hydrogen 
or helium is the hydrogen Lya line. (Trace metals produce 
additional absorption features, though these can be discrim- 
inated based on linewidth or because they have a doublet 
feature, e.g., C iv.) The He I band is the analogous region for 
He I, where the He I 584 A lines are the only source of he- 
lium opacity; however, this region also may contain absorp- 
tion due to all of the Hi Lyman transitions, confusing the 
assignment of any individual line to He I. In principle, unsat- 
urated He I lines could be identified because they are twice 
as narrow as unsaturated Lya lines. However, line blend- 
ing and the expected weakness of the He I lines complicate a 
practical implementation of that criterion, and spectroscopy 
of sufficiently high resolution would require an extremely 
bright target. 

Even after hydrogen reionization is complete at 2 ~ 6, 
the filamentary nature of the clumpy IGM leads to signifi- 
cant Lya opacity in regions with neutral hydrogen column 
density of A^hi ~ 10^'* cm~^. A line of sight through the uni- 
verse pierces many of these filaments at different redshifts, 
resulting in the Lya forest described above. Helium is ex- 
pected to be singly reionized (24.6 eV ionization threshold) 
at a similar redshift to hydrogen, for almost any expected 
ionizing spectrum. Consequently, absorption by neutral he- 
lium at 2 < 6 should be confined to the same filaments that 
give rise to hydrogen absorption.^ This justifies a one-to- 
one search for He I 584 A lines at the same redshift as lines 
identified in the 'Lya band' region, where the Lya band is 
defined as the spectral region where absorbers with a He I 
584 A line in the He I band exhibit Lya absorption (see Fig. 

IB 

The application of absorption line studies of the IGM 
is limited by the supply of suitable background sources, pri- 
marily quasars. The Sloan Digital Sky Survey^ (SDSS) is a 
large photometric and spectroscopic survey of the northern 
sky that will identify and measure redshifts for almost all 
quasars brighter than 20th magnitude over one quarter of 
the sky. At the current discovery rate of the survey, it will 
find approximat ely 1000 quasars at 4 < 25.2, and about 200 
at 4.8 < 2 < 5.2 ('And erson et al.l 2001:'Fan ct al. 2001bY As 
we will demonstrate (Section|7|l, rest-frame extreme ultravi- 
olet (EUV) absorption line study of 2 ~ 5 quasars requires a 
large sample of quasars, and could only be realized with the 
large catalog of 2 ~ 5 quasars that the SDSS will provide. 



3 IONIZING BACKGROUND 

Next we consider the relationship between the ionization 
state of IGM absorbers and the ionizing background. The 
quantity of interest for an absorbing system is the ratio of 
column density of neutral helium, A^hoI, to the column den- 
sity of neutral hydrogen, A'hi, 



(1) 



As we will describe in Section|l] the abundance of absorbers 
as a function of redshift and A^hi is already observed. Thus 
predictions for He I absorption from these systems can be 
easily made from calculations of the relative ionization states 
of helium and hydrogen. 

First we assume that the relative abundance of neutral 
helium to neutral hydrogen in an absorbing system with 
high Hll and Hell fractions'^ is described by the optically- 
thin limit (e.g.. lMiralda-Escude fc OstrikeI^^1992^ ■ 
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where whc/wh — 0.083 is the cosmic abundance ratio of 
helium to hydrogen (assuming a mass fraction of *He of 
0.25, consistent with the combination of WMAP CMB re- 
suhs (SDerEcl et al.'ioO?) and big bang nucleosynthesis pre- 
dictions (.Buries et al. 2001i'l 'l. aHci/ani = 1.0 is the ratio 
of the radiative r ecombination coefficients for He I and H I 
llOsterbrocklll989l) . and 
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(3) 



and similarly for FhcI. The angle- averaged specific intensity 
is Jv, cri^v) is the photoionization cross-section of species i, 
and the lower limit of the frequency integral is i^i, the photon 
frequency at the ionization threshold of species i =H,He. 

The most important absorbers for measuring r) are sys- 
tems where the observed equivalent widths of the Lya forest 
absorption lines are roughly equal to the spectral resolution 
obtained in the Lya forest spectrum. For a typical spectral 
resolution of 0.27 A (a choice justified in Section [6.21 . this 
corresponds to A^hi — 3 x lO'^* cm~^ at 2 ~ 5. At such col- 
umn densities, the optical depth to photons at the hydrogen 
ionization threshold is thi J; 2x 10"'', thus the optically-thin 
approximation should be valid. 

However, optically thick absorbers can have a substan- 
tial impact on the spectrum of the ionizing background, both 
through absorption and emission. Haardt & Madau ( 199^ 
find that the mean UV background is attenuated by roughly 
a factor of 3 at both i^hi and t'Hei; additionally, much of the 
flux shortwa rd of 228 A is absorbed by He 11 and re-emitted 
at 304 A (see lHaardt fc Mada u" 1996'. fig. 5c). These changes 
to the unabsorbed spectrum have important implications for 
the Hem abundance, but do not affect the abundance ratio 
of He I to Hi. 



^ This is in contrast to doubly ionized helium, which is dis- 
tributed more diffusely (that is, the absorption featur es traces 
smaller over-densities) than hydrogen iKriss et al.1l20o'Ill . 
^ http://www.sdss.org 



^ We neglect the possible presence of He III in the absorption sys- 
tems. Observations sugges t that the second reionization of helium 
did not occur until 2 3 <Kriss et alJl200lHTheuns et al.ll2002^ . 
and IWvithe fc Loed <2003lj h predict that the full helium reion- 
ization happens over a narrow redshift interval. These studies 
suggest that most of the intergalactic helium was singly ionized 
at z ~ 5. 
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Absorption by Hi is quantified by discretizing the ab- 
sorption into individual absorbing systems, identified by 
their redshift, z, and neutral hydrogen column density, A^hi- 
The distribution of H I systems may be described reasonably 
well by power laws in both redshift and column density. We 



0.4 
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Figure 2. The ratio, r/thim of the column density of neutral he- 
lium to the column density of neutral hydrogen, as a function 
of the spectral slope of the ionizing background, cn,. We assume 
that the absorbing systems are optically thin to ionizing radia- 
tion, highly ionized in hydrogen and highly singly-ionized in he- 
lium. The effective power-law index of a normal stellar population 
is O]-, = 2.05; the effective power law index of a population of very 
massive metal-free stars is cn, = 0.28. 



Because we are only concerned with the relative rate of 
neutral-helium ionizations to hydrogen ionizations, we ap- 
proximate the ionizing background by a power law. 



(4) 



Using a continuous star-formation history with metallic- 
ity Z = 10"^ and a Salpeter IMF from 1 to 100 Mq 
jLeitherer et alJll999^ . the effective power-law index is ab = 
2.05. Also of interest is the spectrum of very massive metal- 
free stars, which represent the most extreme expectation of 
the stellar IMF at very high redshift: these stars are well 
approximat ed by a blackbody spectrum wit h temperature 
T = 10^ K jBrornm. Kudritzki fc Loebll200lll . Very massive 
metal-free stars yield a particularly hard power-law index of 
Qb = 0.28. Quasars produce a spectrum with a power-law 
shape matching their own EUV spectral slope; that value 
is not well constrained at high redshift, but Qb = 1.6 may 
be a reasonable guess (see Section [SJl . Finally, x-rays r adi- 
ated from early stars have a near-zero spectral slope Johl 
I2OOIII . but the effect of secondary ionizations should pro- 
duce an effective slope slightly greater than zero (S.P. Oh, 
private communication) . Figure |5| shows the value of r^thin 
as a function of Qb. 



4 IGM ABSORPTION MODEL 

Our model for the Hi absorption in the IGM is empirical, 
based on the data obtained in many quasar absorption line 
studies. Ongoing and planned observations are likely to con- 
strain its parameters better in the future. 



systems, and systems with Ahi ^ 1.6 x 10 cm~ as Ly- 
man Limit Systems (LLSs) . The distribution of H I systems 
in redshift and column density are parametrized by three 
constants: A, 7, and s, 



dNjz) 
dz 



A{1- 



/(Am) « Ah 



(5) 
(6) 



Here N{z) is the total number of absorbers along a line of 
sight to redshift z, and /(Ahi) is the number of absorbing 
systems per unit H I column density. 

We chose s = 1.5 for the power law index of 
the Ahi distribution for all values of Ahi, and assume 
that the a bsorber population extends u p to Ahi = 
10^^ cm-2 JStorrie-Lombardi fc Woi^l2000D . For the Lya 
forest, we adopt the redshift evolution det ermined by H ST 
observations at low redshifts ibobrzvcki et all l2002l . cf. 
IWevmann et al .' 1998) and from ground-based surveys at 
high redshifts 7Bechtoldil99il . Those surveys covered the H I 
column density range of lO'^* cm~^ Ahi ^ 1.6 x lO^'^ cm"'^ 
(i.e., from a rest-frame equivalent width of 0.24 A up to the 
LLSs). We extend their results to Ahi = 10^^ cm~^ by scal- 
ing up the number density using the Ahi distribution quoted 
above, to give 



A = 106, 7 = 0.65, for z < 2.4, 
A = 29, 7 = 1.7, for z > 2.4. 



(7) 



The transition redshift, z — 2.4, was chosen to match the 
low-redshift and high-redshift fits smoothl y; this is higher 
than the traditional value of about 1.5 (e.g.. lWevmann et alJ 
I l99a)i but is a conseq uence of the upward revision of 7 b y 
iDobrzvcki et alJ (l2002l) compared to lWevmann et all il99d) ; 
given the large uncertainties in the values of A and 7, this 
is not alarming. 

For the LLSs, we adopt A = 0.2 and 7 = 1.5, which 
are co nsiste nt with the results of IStorrie-Lombardi et alJ 
('1994*) a nd IStengler-Larrea et alJ ^1995^ . We note that 
Stcnglcr- Larrea et^T^995^ do "not present an analysis of 
their full sample of LLSs because they exclude absorption 
systems within 5000 km s~^ of th e quasar, yet their re- 
sults a re consistent with the result of lStorrie-Lombardi et all 
il994f) for all LLSs. These parameter choices provide a rel- 
atively smooth intersection between the abundance of LLSs 
and LyQ forest systems at 2 ^ 1, assuming the column 
density distribution described above. The covariance of the 
uncertainties in A and 7 implies that the mean number 
of LLSs toward a high-redshift quasar is constrained, but 
the LLS distribution with red shift is uncertain: the 1-a un- 
certainty on 7 is about 0. 4 iStorrie-Lombardi et al.l[l994l : 
IStengler-Larrea et alJll995ll . 

The absorber distributions described above can be con- 
verted into a mean comoving H l number density using 



(n^i(z)) = A(l + z)- 



HoEjz) 
c{l + zy 



{Nm) , 



(8) 
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where Hq is the Hubble constant, (A^hi) is the mean Hi 
column density of the absorber population, and 

E{z) = + z)^ + s^A + (1 - n,n - nA){i + zf] '''^ . (9) 

For a universe with density parameters f2m = 0.3 in 
matter, Qa = 0.7 in a cosmological constant, and 
Qhh^ = 0.02 in baryons, and a Hubble constant h = 
//o/(100 km s-^Mpc"^) = 0.7, this gives an H I neutral frac- 
tion at 2 > 4 of, 

(xHiW) ^ 9.8 X 10-^A(1 + .)-V2_J^gil_. (10) 

For LLSs (which have {JVhi} = 4 x 10^^ cm'^), A = 0.2 
and 7 — 1.5, thus (xiii( z = 6)) = 0.01. T his is consistent 
with the limit derived bv lFan et al.l i2002f) from Hi absorp- 
tion toward a, z = 6.28 quasar. 

We assume that there is no correlation between absorb- 
ing systems. There is obse rvational evidence for clustering of 
the L yg forest lines fe.g.. lLiske et al.ll200ll : IPobrzvcki et alJ 
|2002| . and references therein), but our assumption is con- 
servative for the purposes of this paper: we will show that 
to study He I absorption, we need quasars with fewer than 
average absorbers along the line of sight. Clustering of ab- 
sorbers would skew the 'bad' quasar targets (i.e., those with 
many strong absorbers) worse, and skew the 'good' quasar 
targets better. This is not an important effect for Iow-A^'hi 
absorbers, since they are very numerous, but could be a sig- 
nificant effect for LLSs if they arc clustered (note, however, 
that Sargent, Stcidcl, & Bokscnbcre 1989 showed that LLSs 
followed Poisson statistics in their sample of 37 absorbers). 
The proximity effect generates an absorber deficit of ~ 30 
per cent within 4 Mpc of the quasar, or about 2 00 A 
observed from the quasar Lya line JScott et al.ll2000D . We 
also ignored this relatively small effect. 

In addition to the absorber redshift and column-density 
distributions, our IGM model must also describe their ab- 
sorption proper ties. We used the photoioni zation cross- 
sectio ns given bv lOsterbrock * (1989) for H I andl Verner et all 
il99dh for He I. Our model treats the first 11 line transitions 
from the ground state (which all atoms are assumed to oc- 
cupy), using data from the NIST Atomic Spectra Database 
V2.0.'' The line profiles were modelled as Doppler cores of 
width 6 = 26 km s~^ \Kim e t al. 1997) with damping wings 
outside of the core (e.g.. IPeeblea, 1993.) .^ 

Our absorption model includes only hydrogen and he- 
lium; metal lines are not considered here. There are not likely 
to be many metal line systems in the quasars used to mea- 
sure r] because metal lines are associated with H I absorbers 
of high column density, and, due to the selection techniques 
employed (see Section|SJ , quasars with strong absorbing sys- 
tems are unsuitable for measuring rj. Moreover, the strongest 
metal lines that could pollute the Lya band, such as the C iv 
and Mgll doublets, show only a few absorbers per unit red- 
shift, compared to of order 100 Lya absorbers in the same 
redshift interval. In the He I band, metal lines are again far 
less numerous than H I lines, and thus were ignored. 

* http:/ /physics. nist.gov/cgi-bin/AtData/main_asd 
^ In Section l6.2l we assume instrumental resolutions considerably 
larger than 26 km s~^, so the choice of exactly that value and 
our use of a single b value for all absorbers rather than the actual 
distribution should have a negligible effect on our results. 



One complication to the application of absorber statis- 
tics determined from other quasar samples to the population 
of quasars discovered by SDSS is the quasar selection tech- 
niques employed. The SDSS quasar selection primarily uses 
observed optical colors to generate a well-defined selection 
function jFan et alJIlioil^ : by contrast, the quasars studied 
in absorption line surveys we re culled from heterogeneous 
catalogs (e.g.. .Hewitt fc Burbidge 19871 . and include radio- 
, x-ray-, and emission-line-selected quasars in addition to 
color-selected quasars. An analysis of the systematic errors 
introduced by using absorption lines in quasars discovered 
by several different techniques to predict absorption patterns 
in quasars selected by another technique is outside the scope 
of this paper; we simply analyze the colors of the quasars 
most important to this study to ensure they would meet the 
quasar color-selection criteria of SDSS (see Section |7|l. 



5 MODEL QUASAR INTRINSIC SPECTRUM 

In order to simulate realistic observations of IGM absorption 
of background quasars, we need to make assumptions about 
the intrinsic properties of the quasar UV spectrum. The 
SDSS will provide a large catalog of high-redshift quasars, 
so we attempt to model the typical properties expected of 
SDSS quasars based on the results of the survey so far. 

In a sample d rawn from 182 square degrees of the SDSS, 
iFan et al] (|20^^ presented 18 quasars with z ^ 4 and i* < 
20, where i* is the preliminary SDSS determination of the 
quasar i' magnitude. We assume the quasars presented are 
a fair sample of the final SDSS resu lts. ^^^^ 

For each quasar in their sample. IPan et aD ll2001all com- 
pute the properties of the quasar near-UV (NUV) contin- 
uum, assuming a power law form, oc ;/~"nuv They set 
the normalization to the quasar continuum at the observed 
wavelength corresponding to rest-frame 1450 A. They then 
convert that specific flux into a magnitude on the AB sys- 
tem, 

ABi450 = -2.51ogio[/4A,eBt = 1450 A)] - 48.6, (11) 

where /,y(Arcst = 1450) is the specific flux in 
erg s"'^ cm~^ Hz~^ at rest-frame 1450 A. Typically 
ABi45o ~ i* or z* , depending on quasar redshift. iFan et all 
t2001ah also estimate qnuv for each quasar; the mean for 
quasars with 2: ^ 4 is 0.6, with substantial u ncertainty in 
the slope of any individual quasar. However, iTelfer et al.l 
(2002) found a break in the UV slope of their quasar com- 
posite spectrum near the wavelength of Lya; the measured 
slope be came softer blue- war d of L ya. This confirmed earlier 
work bylzhen^eTaO The iTelfer et all (120021) sam- 

ple is comprised almost entirely of quasars with z < 2.5; 
for radio-quiet quasars they find a mean EUV slope of 
about 1 .6, with a break to a shallower NUV slope at about 
1250 A. ITelfer et all (l2002f) compare thei r radio-quiet sample 
to the radio-quiet SDSS quasar sample of lVanden Berk et alJ 
i200lll . and find an evolution toward harder NUV slope with 
samples at higher redshift; EUV comparison between the 
samples is extremely difficult due to the presence of strong 
IGM absorption in the high-redshift quasars. 

Given the available observational evidence, we adopt a 
conservative model for the EUV properties of SDSS quasars: 
AB1450 = 20 with a NUV slope of 0.6 between rest-frame 
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1450 and 1250 A and an EUV slope shortward of 1250 A 
given by oeuv = 1-6. 



6 LINES OF SIGHT TOWARD SDSS QUASARS 

Many z ~ 5 quasars will be discovered by SDSS; in this 
Section we evaluate how good the best quasar target will 
be. We then simulate observations of the best quasar target 
and demonstrate our ability to constrain the spectrum of 
the ionizing background from such data. 

Extrapolating from the results of the color-selected sam- 
ple of Fan ct al. ( 2001a), SDSS wiU discover ~ 1000 quasars 
with 4 < z < 5.2. We assume the SDSS quasar selection 
function does not depend much on redshift over the range 
4 < z < 5.2; this is a conservative estimate for our pur- 
poses because the detection probability is higher f or quasars 
with 4 . 7 < z < 5.2 than f or other z > 4 quasars iFan et alJ 
|200li).|^nerai] ^2001b^ fit the redshift dependence of the 
high redshift quasar spatial density with p oc 10~°'^^. Thus 
about 200 of the high-redshift quasars will fall within the 
range 4.8 < z < 5.2. 

We would like to know how many of these quasars have 
spectra like the middle panel of Fig. which would be 
useful for measuring He I lines, and how many have spec- 
tra like the bottom panel of Fig. which would not be a 
suitable target for measuring He I lines. We quantify this 
by determining the distribution of quasar sightline 'suit- 
ability,' as measured by the flux of the quasar in the He I 
band. From that distribution and the expected number 
of SDSS quasars, we determined what the most suitable 
quasar discovered by the SDSS will be. Though for some 
applications an analyt ic approach to the subject is suitable 
JZuo fc Phinnevll9^3^■ we require a Monte Carlo simulation 
approach jMoller fc JakobsenlllOQol : IJakobsenlHool) for this 
study. 

6.1 Absorber Monte Carlo simulations 

We started from the distribution of H I absorbing systems 
in redshift and column density described in Section |1| For a 
quasar at a given redshift zq, we used Poisson statistics to 
generate the number of low-z Lya forest, high-z Lya forest 
and LLS absorbers along the line of sight. Once we generated 
the number of absorbers of each type, we assigned redshifts 
and column densities drawn from the absorber model de- 
scribed in Section |1] This list of absorbers was then passed 
through a routine to calculate the optical depth to H I at ev- 
ery sampled wavelength value. ^ In each simulated spectrum, 
the optical depth due to helium was stored separately from 
hydrogen; given the universal optically thin ratio of column 
densities, r^thin derived in Section|3] the helium optical depth 
scales simply as rythin. 

For each LOS, the Hi optical depth data were con- 
verted into the IGM transmission as a function of wave- 
length. We integrated the IGM transmission over the He I 

^ Wavelength sampling ranges from 0.08 to 0.28 A, chosen so 
that there are three samples per Doppler FWHM of the spectral 
lines of interest (He I lines are twice as narrow as Hi lines); all 
results were checked for sensitivity to spectral resolution and are 
converged at this resolution choice. 
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Figure 3. Cumulative histogram of the Hel-band transmission 
along quasar lines-of-sight, as a function of redshift and absorber 
model. The solid curve shows the fraction of z = 5 quasars that 
have a mean transmission in the He I band greater than a given 
mean transmission value T(Hel), assuming s = 1.5 in the dis- 
tribution of absorber column densities (Section The lowest 
dashed curves are for s = 1.4 and s = 1.6, respectively, assum- 
ing 2 = 5 quasars. The two dotted curves closely above the solid 
curve are for z = 4.9, 4.8, and the two dotted curves closely below 
the solid curve are for z = 5.1, 5.2, all assuming s = 1.5. The hor- 
izontal dashed line at 1/200 intersects the curves at the expected 
T(Hel) of the best quasar in a sample of 200 quasars. 

band to derive the mean Hel-band transmission along each 
LOS, r(Hel). Figure shows the cumulative fraction of 
quasars with T(Hel) greater than a given value. For sim- 
plicity, we will assume the expected 200 SDSS quasars with 
4.8 < ZQ < 5.2 aU faU at zq = 5. Then from Fig. El we 
can read off the value of T(Hel) corresponding to a frac- 
tion of 1/200; this is the largest expected r(Hel) value of 
the lines-of-sight toward the SDSS quasars. The result is 
T(Hel) — 0.03. Note that if we change the assumed red- 
shift of the quasars by Az — 0.2 (the dotted lines in Fig. 
1^, the expected maximum value of T(Hel) for 200 quasars 
changes by a factor of about 2. When the final SDSS quasar 
catalog is constructed, we may use the real distribution of 
quasar redshifts together with our simulation machinery to 
generate a r(He l) histogram that depends only on the IGM 
absorber model. Thus measurements of the observed distri- 
bution of T(Hel) can be used to constrain a combination 
of the properties of the H I component of the IGM and the 
intrinsic quasar spectral shape. 

Small changes to the absorber-model parameter s, the 
slope of the column-density distribution, have a strong effect 
on the expected maximum v alue of r(He l) (as pointed out 
bv lMoUer fc Jakobsen|[T990l . and illustrated by the dashed 
lines in Fig.jSJ. In the next subsections we discuss what could 
be learned about the z ~ 5 ionizing background under the 
assumption that SDSS will discover a quasar with properties 
given in Section |K| along a LOS with T(Hel) = 0.03, our 
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'expected' value for the best SDSS quasar, or r(Hel) — 0.1, 
which we consider to be a reasonable 'optimistic' value for 
the best SDSS quasar. 

6.2 Quasar absorption line spectrum simulations 

Based on the results of the previous subsection, we ana- 
lyzed lines-of-sight with r(Hel) — 0.03 and 0.1, our 'ex- 
pected' and 'optimistic' values for the best SDSS quasar 
LOS. For each observational realization of a quasar along 
the simulated LOS, a value of rj was assumed. Then the 
optical depth data were used to construct the IGM trans- 
mission as a function of wavelength. This was multiplied by 
the intrinsic quasar spectrum from the model described in 
Section]^ to generate the obser ved spectrum. The Galacti c 
extinction curve was taken from lSchefBer fc ElsasseJ il98d) . 
Atmospheric extinction was applied based on the sky trans- 
parency from Mauna Kea (the location of the Keck Ob- 
servatory).'^ Our modelling of the telescope and instrument 
parameters included slit losses, system throughput, instru- 
mental resolution, pixel sampling, spatial extraction of the 
spectra, and detector noise. Slit losses and spatial extraction 
both incorporate an assumed value for the seeing. A mean 
emission spectrum of the sky,® measured at Mauna Kea, was 
also passed through the telescope/instrument model. 

Our model observations were performed assuming a 
Galactic extinction value of A„/ = 0.15 toward the t arget 
quasar, typical for SDSS quasars JSchneider et alJ2002ft . We 
assumed the observations were made at a constant value of 

I. 2 airmasses and a constant seeing of 0.7 arcsec. 

Our simulated observations of the Hi Lya forest were 
based on the properties of the Keck II 10-meter tele- 
scope and the Ec hellette Spectrograph and Imager (ESI; 
ISheinis et al.ll20o3) , an intermediate resolution optical spec- 
trograph, with ESI in echelle mode (R ~ 4000). We as- 
sumed a 1 arcsec slit, with 75 km s~^ resolution and 

II. 5 km s~^pixel~^. The spatial extraction window was 
0.77 arcsec; the spatial pixel scale was 0.154 arcsec pixel"^. 
We assumed no dark current, and a readnoise of 2.1 counts 
per pixel (the gain was 1.3 e~ per count). We assumed each 
Hi Lya forest spectrum consisted of 40 co-added 1000 sec 
observations, for a total exposure time of 11 hours. This is 
similar to the exposure times of Keck observatio ns of the 
highest-redshift SDSS quasars (IWhite et alj|2003ll . 

We assumed that the He I band was observed with the 
Keck I 10-meter tele scope and the Lo w-Resolution Imaging 
Spectrometer fLRIS: IOke et al.lll995H using the blue chan- 
nel of LRIS fLRIS-E nMcCarthv et allll998l) with the 1200 
line grism. This configuration with a 0.7 arcsec slit deliv- 
ers a resolution of 1.30 A, sampled at 0.24 Apixel"^. The 
spatial extraction window was 1.08 arcsec, about 1.5 times 
the seeing; the spatial pixel scale was 0.2151 arcsec pixel"^. 
We assumed no dark current, and a readnoise of 2.5 counts 
per pixel (the gain was 1.6 e~ per count). Because readnoise 
otherwise made a non-negligible contribution to the noise, 
we binned the data by 2 pixels in the spatial and spectral 

http:/ /www2. keck. hawaii.edu/inst/lris/atm_trans. html 
* http:/ /www. cfht.hawaii.edu/Instruments/ObserverManual/ 
chapter5.html; note the sky will darken slightly as we approach 
solar minimum in 2008. 
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Figure 4. A simulated absorption spectrum for a 2 = 5 quasar 
with AB1450 = 20. The top panel zooms in on the 2 = 4.8-5 Hi 
Lya forest of the spectrum from the middle panel of Fig. The 
second panel shows the corresponding 2 = 4.8—5 He I 584 A forest 
region of the same spectrum; the dotted curve is the absorption 
from Hi only, and the solid curve assumes Jjthin = 0-5 for illustra- 
tion. The lower two panels show mock observations of the quasar, 
using the observation model described in Section 16.21 

directions. We assumed each He I 584 A forest spectrum con- 
sisted of 72 co-added 2000 sec observations, for a total expo- 
sure time of 40 hours. This is substantially more observation 
time than is typically devoted to any one object with Keck. 
However, the total observation time devoted to study Hell 
absorption in the z — 2.885 quasar HE2347-4342 was almost 
150 hours with the FUSE satellite. 

For each simulated exposure, we computed the expected 
signal, expected sky, and expected detector noise. The ob- 
served counts in each spectral pixel due to source, sky and 
detector were drawn from Poisson distributions with the 
means set to the expected counts, and the three contribu- 
tions were summed. Then the expected sky and detector 
noise were subtracted (assuming photon-noise limited sky 
subtraction). The simulated exposures were then summed 
to create the final simulated observation. 



6.3 Quasar absorption line spectrum analysis 

Our goal is to estimate the ionization state of the IGM as 
measured by the relative abundance of He I to Hi, rj. Our 
analysis technique is to make use of the cross-correlation of 
the H I Lya forest spectrum and the He l-band spectrum, 
which depends on rj as illustrated in Fig. 2] 

Each pixel in the Hel-band spectrum was matched up 
with the Lya-band spectrum pixel whose wavelength is clos- 
est to Am = (1216/584)AhoI, where AhoI is the wavelength 
of the He i-band pixel. Thus an absorber at a given redshift 
will exhibit He I 584 A absorption and Lya absorption in 
matched pixels. Figure 13 shows a plot of the fluxes, normal- 
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Figure 5. An example of the pixel-by-pixel cross-correlation of 
the absorption in the He I band and the H I Lyo band. The solid 
squares represent continuum-normalized flux pairs from an obser- 
vation of a quasar LOS with r) = 0.06; the open squares are from 
an observation of the same LOS but assuming ri = 0.20. 



ized to the local quasar continuum, of pixel pairs matched 
in this manner, for two values of r]. If He I 584 A absorption 
and LyQ absorption were the only features in the spectra, 
and the observations were perfect, then the points would 
fall along a tight locus described primarily by rj. In regions 
of little absorption, both continuum-normalized pixel fluxes 
would be near one: regions of strong absorption would have 
normalized Hi pixel fluxes near zero and normalized He I 
pixel fluxes set by -q. However, the presence of other lines 
and the effects of instrumental smoothing and observational 
noise scatter the points. There is little scatter in the Hl- 
pixel-flux direction, both because these simulated observa- 
tions have a good signal-to-noise ratio, and because there 
are no absorption features in the Lya band besides LyQ. 
There is substantial scatter in the He l-pixel-flux direction, 
though, because of lower signal-to-noise in the observations, 
and also the presence of Lya absorption (and, to a lesser ex- 
tent, absorption due to the other Lyman series lines) from 
low-redshift absorbers. Despite the scatter, one can easily 
see that, at low values for the normalized Hi pixel flux, the 
points Fig. |5] for the r) — 0.06 case (solid squares) lie higher 
(less He I absorption) than the points for the i] = 0.20 case 
(open squares). 

In the He l-band spectrum, we would like to know what 
the quasar effective continuum level is after accounting for 
bound-free absorption by higher-redshift absorbing systems. 
For a zq = 5 quasar, any strong absorber between z = 2.53 
and z = 2.84 will cause a change in the effective continuum 
in the He I band, due to the Hi bound-free absorption edge 
at 912(1 -I- z) A. Consequently, the Hi Lya forest spectrum 
should be searched for strong Lya lines over the range of 
4292 to 4669 A (for = 5). We assumed for our analysis 
that the effective continuum can be accurately estimated. In 



Figure 6. Contour plot of the ensemble continuum-normalized 
fluxes of pairs of corresponding He l-band and Lya-band pixels 
(see Section 16.31 . The solid contours assume rj = 0.06, and the 
dotted contours are for rj = 0.20. Pixel pairs with strong Hi ab- 
sorption have H I pixel fluxes near zero; for those pairs increasing 
r] shifts the He: pixels to smaller fluxes, due to increased He I 
absorption associated with the H I absorption. 



the analysis of real observations, one could apply the same 
technique used to estimate the effective continuum to our 
mock observations, and thus analyze the real observations 
in an unbiased way. 

Using our large Monte Carlo library of simulated quasar 
LOSs, we generated the distribution of pixel pair fluxes, as 
a function of rj, for multiple observational realizations of a 
quasar along each LOS. Figure |S| shows contours that rep- 
resent the likelihood of finding a given pair of continuum- 
normalized fluxes for 77 = 0.06 and r; = 0.20. This distribu- 
tion of mock-observational results were then used to analyze 
an individual observation of a quasar LOS: we generated 
the matched pair data as shown in Fig. |S] then summed the 
value of the likelihood at that point over all the pixel pairs 
with continuum-normalized Hi pixel fiuxes less than 0.5 (cf. 
Fig-injj for each of the two values of 77 under consideration. 
We termed each sum of the likelihood values (which depend 
on ri) Lijj). The L{rj) value is larger when a set of pairs 
matches up well with the contours of Fig.|S]for 77, and smaller 
when there is not a good match. Then we formed the ratio 
R = L(0.06)/L(0.20), which we find is better for recovering 
our input value of r) than, for example, L(0.06) — L(0.20). 
We have found that binning the He l-band spectra by 2 pix- 
els (in addition to the 2x2 on-chip binning) reduces the 
scatter in R. 

FigureQshows histograms of R for each of our two input 
values of rj, assuming observations described in Section r6.2l of 
a quasar along a sightline with r(Hel) = 0.03, our expected 
value for the best SDSS LOS toward a z ~ 5 quasar (see 
Section [6. 11 . The histograms are clearly separated, though 
they do show some overlap. Figure |H| is the same as Fig. |7| 
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Figure 7. Histograms of R (see Section l6.3i for each of our two 
input values of -q, 0.06 and 0.20, assuming observation of a quasar 
along a sightline with T(Hel) = 0.03. The solid (dotted) line 
shows the R histogram for analysis of quasar LOSs with rj = 0.06 
(r; = 0.20). 
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Figure 9. Normalized cumulative histograms of R from Fig.|7| 95 
per cent of the LOSs have i?(0.06) > 0.958, the position of the left 
vertical line, and 95 per cent of the LOSs have i?(0.20) < 0.993, 
the position of the right vertical line. Thus if the intrinsic rj value 
along a LOS is either 0.06 or 0.20, one observation has about a 
60 per cent chance of rejecting the other value. 
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Figure 8. Same as Fig.[3but for an LOS with r(Hel) = 0.1. 



but for observations along a LOS with r(Hel) — 0.1, our 
optimistic value for the best quasar LOS in the SDSS. In 
this case the histograms are much more clearly separated, 
illustrating the importance of discovering a very good quasar 
LOS. 

Figures |U] and IIUI cumulate (and normalize) the his- 
tograms of Figs.[7|and|H| They illustrate that, for our model 
observations toward one r(Hel) = 0.03 LOS, if the true 



value of Tj is 0.06, then we have almost a 60 per cent chance 
of making an observation that will reject rj = 0.20 at the 
95 per cent confidence level. Conversely, if the true value 
of ri is 0.20, then we have a 60 per cent chance of making 
an observation that will reject ri = 0.06 at the 95 per cent 
confidence level. The situation along a r(Hel) =0.1 LOS is 
much more optimistic: if the true value of rj is either one of 
our model choices, we will reject the other value of rj about 
95 per cent of the time at the 95 per cent confidence level. 

Using the observations we propose, it will be difficult 
to estimate a precise value of from the data. If the true 
value of rj is, for example, 0.11, then we will only be able 
to reject very extreme values of rj. However, if the ionizing 
background is either hard or soft, it may well be possible 
to reject the other hypothesis. For example, if we measured 
an R value of 1, we could be confident that the ionizing 
background at z ~ 5 is not dominated by Pop ll stars. 



7 OBSERVATIONAL QUASAR SELECTION 
TECHNIQUES 

In previous Sections we estimated the likely properties of the 
best SDSS quasar for the measurement of He I absorption 
features, and what we could learn about the z ~ 5 ionizing 
background from analysis of such a quasar. In this Section 
we describe the final aspect of practical implementation, how 
to pick out the best SDSS quasar from the expected sample 
of 200. 

We would like to select the 2 ~ 5 SDSS quasar with 
the highest flux in the He I band. So far we have primar- 
ily discussed T(Hel), the IGM-transmitted fraction in the 
He I band; however, this was always under the assumption 
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Figure 10. Normalized cumulative histograms of R from Fig. ISl 
95 per cent of the LOSs have _R(0.06) > 0.970, the position of 
the left vertical line, and 95 per cent of the LOSs have ij(0.20) < 
0.972, the position of the right vertical line. Thus if the intrinsic 
rj value along a LOS is either 0.06 or 0.20, one observation has 
about a 95 per cent chance of rejecting the other value. 



of the intrinsic quasar model described in Section |3 The 
real SDSS quasar sample will certainly include quasars with 
a range of values for AB1450, and these quasars may have 
a range of EUV spectral shapes. Due to these distributions, 
the quasar with the largest Hel-band flux may not be the 
quasar with the largest r(Hel); we made conservative as- 
sumptions about AB1450 and the EUV spectrum, so we ex- 
pect that the true distributions of quasar properties may 
only improve the suitability of the best SDSS quasar over 
our estimates. 

For ZQ — 5, the He I band covers 3222 to 3506 A; this 
is withi n the SDSS u' fi lter, which runs from ~ 3250 to 
3750 A llFan et alJl2001ah . As a consequence we expect the 
u' magnitude to serve as a good proxy for the Hel-band 
flux. Figure [TTI shows the relationship between r(Hel) and 
T(m'), the transmitted flux fraction in the u' band. The 
strong correlation illustrated in the scatter plot is indepen- 
dent of the intrinsic quasar properties. Assuming our stan- 
dard quasar model, the transmitted fractions are directly 
related to fluxes; the cumulative histogram at the bottom 
of Fig. 111! shows the fraction of observed quasars brighter 
than a given value of u' . We expect that 95 per cent of 
SDSS « ~ 5 quasars will be fainter than u — 27, but the 
quasar with the highest Hel-band flux will have u' ~ 25.5, 
with an optimistic value of u = 24. The u'-band l-cr limit- 
ing magnitude for the SDSS survey is expected to be about 
u' — 24 (.Fan et al. 2001a) . thus under optimistic assump- 
tions there is some chance the SDSS survey itself would make 
a u detection of a 2; ~ 5 quasar, and thus identify an excel- 
lent candidate for spectroscopic follow-up observations. It 
is likely, however, that deeper u'-band photometry on the 
SDSS quasars will be desirable to locate the few brightest 



Figure 11. Plot of T(Hel) versus T{u'). This plot shows the 
mean transmissions in the He I band and in the u' band for 
each z = 5 quasar LOS in 10^ LOS Monte Carlo simulations. 
The excellent correlation between T(Hel) and T{u') at the high 
T(Hel) end of the distribution shows that n'-band photometry 
is an effective method for selecting the best quasar for follow-up 
He l-band spectroscopy. Under the assumption that the intrinsic 
quasar spectrum is given by the model of Section |5] we convert 
T(ji') into a -u'-band AB magnitude, and plot the cumulative 
fraction of quasars brighter than that magnitude in the bottom 
panel. The horizontal dashed line at 1/200 intersects the curve at 
u' = 25.3, the expected u' magnitude of the brightest quasar in 
a 200 quasar sample. 

quasars in the u band. The 2.5-meter SDSS telescope inte- 
grates for only 54.1 seconds per filter, so achieving deeper 
u photometry on a small telescope equipped with a blue- 
sensitive CCD should not be difficult (extending the l-cr 
limit to u' = 25.5 would only require 2.5 minutes per source 
with the SDSS telescope, or about 8 hours of total integra- 
tion time for all 200 SDSS 2; ~ 5 quasars). 

Once candidates have been identified from it'-band pho- 
tometry, low-resolution spectroscopy should be performed 
to determine the flux in the He I band; because the u' band 
extends longward of the He I band, a quasar can be rela- 
tively bright in u but faint in the He I band (note the points 
scattered below the main correlation in Fig. IIH . The low- 
resolution spectroscopy will require a blue-sensitive spectro- 
graph, but on only a 2-meter class telescope. 

We now turn to our final point: the SDSS quasars are 
selected based on photometry in all of the SDSS bands. 
The primary SDSS photometric-selection criteria for 2^ 5 
quasars require u > 22.00 and g' > 22.60 jFan et al.l2001;Jl . 
As illustrated in Fig. 1111 we do not expect any z ~ 5 quasars 
to be nearly as bright as u' = 22.00. In contrast, Fig. 1121 
shows that quasars with r(Hel) = 0.1 will sometimes be 
slightly brighter than 22.60 in g' . Though unlikely, it is pos- 
sible that the SDSS quasar selection algorithm could miss 
not only a 2; ~ 5 quasar, but the best one in the survey for 
our purposes. When the SDSS quasar survey is complete. 
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Figure 12. Plot of T(Hel) versus g' of quasars. These points are 
for the standard intrinsic quasar flux model of Section |5] Only 
quasars with g' > 22.60 (plotted to the left of the dotted line) 
are identified by the SDSS quasar color selection algorithm as 
quasars. Note that this plot focuses on quasars with very op- 
timistic values of T(Hel) > 0.1; quasars with lower values of 
T(Hel) generally all fall at g' > 22.60. 

and the SDSS data are publicly released, a full search of 
the SDSS catalog can be made to determine whether any 
objects meet all the SDSS z ~ 5 quasar criteria except for 
g' > 22.60, but still have, e.g., g' > 22.40. These objects 
could be followed up independently. 



8 SUMMARY 

We have developed a method to infer the dominant stellar 
population responsible for the ionizing background radiation 
at 2 ~ 5, using absorption lines in the spectrum of a z ~ 5 
quasar. Specifically, we related the strength of relative ab- 
sorption in He I lines compared to H I lines to whether Pop ll 
or Pop III stars dominated the ionizing background shortly 
after reionization completed. 

First we related the ionizing background shape to the 
ratio, ri, of He I to H I in an optically thin absorbing system; 
the spectrum of Popll stars gives rise to r; = 0.20 while the 
spectrum of very massive Pop ill stars yields rj = 0.06. Then 
we adopted a model for the population of IGM absorbers 
based on absorption line observations of quasars. The ex- 
isting observations do not tightly constrain absorber popu- 
lations at z > 4; our model does extrapolate to a neutral 
fraction of 1.4 x 10"^ at z = 6.28, consistent with mea- 
surements. We adopted a model quasar spectrum to shine 
through our model IGM; we assumed the quasar flux to be 
at the faint limit for the SDSS quasar survey (20th magni- 
tude in the AB system at rest-frame 1450 A), and modelled 
the EUV spectral shape based on the composite spectrum of 
lower-redshift quasars with a spectral slope of oeuv ~ 1-6. 

We extrapolated the early results of the SDSS quasar 
survey to find that SDSS will discover 200 quasars with 
4.8 < z < 5.2. We used Monte Carlo simulations of our 
IGM model to generate random realizations of lines-of-sight 
(LOSs) toward zq ~ 5 quasars. Based on those results, we 
concluded that a sample of 200 quasars should contain one 
quasar along a line-of-sight with 3 per cent transmission by 
the IGM in the Hel band (537(1 -I- to 584(1 + 2) A); given 
the uncertainty in the absorber model, we consider a LOS 
with 10 per cent transmission to be an optimistic possibility. 

We simulated large numbers of quasar LOSs with 3 and 
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10 per cent transmission. We placed a model quasar at z = 5, 
and then simulated observations of the He I and H I forests. 
Based on the result for a large ensemble of observations, we 
derived an estimation technique that can discriminate (at 
95 per cent confidence) a Pop 11 ionizing background from 
a Pop III ionizing background ~ 50 per cent of the time for 
observations along our expected best LOS, and that can dis- 
criminate the backgrounds 95 per cent of the time for obser- 
vations along our optimistic best LOS. These observations 
require very long integrations with a 10-meter class telescope 
instrumented with a blue-sensitive low-dispersion spectro- 
graph and an optical intermediate-resolution spectrograph. 

Finally, we examined the selection of the best 
quasar/LOS combination on which to perform our proposed 
experiment from the expected 200 SDSS quasars at 2: ~ 5. 
The Hel forest at z = 5 falls in the u' filter; consequently 
u' imaging is an efficient way to select the best quasar/LOS 
combination. However, the SDSS u' photometry will not be 
deep enough to detect 2 ~ 5 quasars (any detection would 
be beyond the optimistic expectation of our IGM model), 
so follow-up u' photometry of the 200 2 ~ 5 quasars will 
be required; these observations will be inexpensive. We also 
consider the possibility that a quasar/LOS combination will 
be so good for our purposes that it falls outside of the 
SDSS quasar color selection criteria. This is very unlikely for 
quasars with ~ 10 per cent or lower transmission, but about 
half of all quasars with 20 per cent transmission, should any 
exist, would be brighter in g' than the SDSS quasar color 
selection limit of 22.6. If quasars are discovered with 10 per 
cent transmission, it may be worthwhile to modify the SDSS 
quasar color selection criteria so as to discover quasars with 
even less IGM absorption in the Hel forest. 
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